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Magneto-optical imaging and Hall-probe array magnetometry are used to measure the field of 
first flux entry, H p , into the same Bi2Sr2CaCu2 08+<5 single crystal cut to different crystal thickness- 
to- width ratios (d/w), and for two angles a between the edges and the principal in-plane crystalline 
(a, b) axes. At all temperatures, the variation with aspect ratio of H p is qualitatively well described 
by calculations for the so-called geometric barrier [E.H. Brandt, Phys. Rev. B 60, 11939 (1999)]. 
However, the magnitude of H p is strongly enhanced due to the square shape of the crystal. In the 
intermediate temperature regime (T < 50 K) in which the Bean-Livingston barrier limits vortex 
entry, there is some evidence for a tiny crystal- orient at ion dependent enhancement when the sample 
edges are at an angle of 45° with respect to the crystalline axes, rather than parallel to them. 

PACS numbers: 74.20.rp, 74.25.Bt, 74.25.Op, 74.25.Qt 



I. INTRODUCTION 



Surface barriers are well-known to delay vortex pen- 
etration into high-T c cuprates (HTSC)i — As a con- 
sequence, the field of first flux penetration H p may be 
significantly higher than the lower critical field H c \ = 
<l>o/47r/ioA 2 In ft (A is the London penetration depth, n = 
A/£ with £ the coherence length, $o = h/2e is the 
flux quantum, and /io = 47r x 10 -7 Hm" 1 ). Among 
possible types of barriers, the Bean-Livingston (BL) 
barrier—— i 22 i 23 arises from the competition between the 
attraction of an entering vortex to the external sur- 
face (this can be described with an image vortex near 
the surface) and its repulsion by the screening (Meiss- 
ner) current. The maximum value of the total free en- 
ergy per unit length due to the introduction of a vor- 
tex, eo In {H c / y/2H a ) , is attained when the latter is situ- 
ated at a distance ^H c /^/2H a from the superconductor- 
vacuum interface. Here, H a is the applied magnetic field, 
H c = \/2ft£o / $o is the thermodynamic critical field, and 
£q = Qq/Att/ioX 2 is the vortex line energy. Therefore, 
for fields close to the penetration field H p = V2H C , the 
BL barrier can be overcome by thermal activation. In 
long thin isotropic superconductors in parallel fields, this 
is thought to happen by the nucleation of a vortex loop 
in the sample bulk££ In layered superconductors such 
as Bi2Sr2CaCu208+5, the vortex loop may encompass 
a single 7 - or several pancake vortices £ Experiments sup- 
porting an important role of the BL barrier in HTSC are 
the observation of asymmetric hysteresis loops accompa- 
nied by magnetic relaxation over the barrier, near the 
critical temperature T c in YBa2Cu307_<5^ and at tem- 
peratures below 50 K in layered Tl2Ba2CaCu208+#^ 
and Bi2Sr2CaCu208+(5.— ^ Further evidence is the ob- 
servation of the predicted sensitivity of the barrier to 
electron^ and heavy-ion irradiation. 8,9 Namely, the pres- 
ence of defects near the specimen surface perturb the 
screening current on the relevant length scale, larger than 



the coherence length £ but smaller than A, thereby favor- 
ing vortex nucleation. 

Following theoretical predictions 24 of a possible 
crystal-orientation dependence of H p in HTSC, the BL 
barrier has recently returned to the center of attention^ 
This orientational dependence would be due to the pre- 
dominant d x 2_ y 2 -symmetry of the superconducting gap 
in HTSC (for a review, see Ref. 1261 ). In that case, quasi- 
particle scattering at surfaces oriented perpendicularly 
to the nodal directions give rise to Andreev bound states 
(ABS)i2£ _ — The quasi-particle current carried by these 
bound states is responsible for the modification of the 
field of first flux penetration.— In Bi2Sr2CaCu208+5 the 
relevant surfaces are those that have a normal oriented 
near to a = 45° with respect to the principal (a, b) crys- 
talline axes^ 

It should be noted that the BL barrier is effective only 
in the situation where the applied magnetic field is par- 
allel to a large flat surface. In the opposite case of a 
superconducting platelet of rectangular cross-section in 
perpendicular magnetic field, flux entry is hindered by a 
barrier of geometrical origin^— The geometrical bar- 
rier comes from the competition between the line tension 
of a vortex cutting through the upper and lower ridges at 
a sample edge, and the Lorentz force due to the Meiss- 
ner current. The penetration field H p is reached when 
the total field at the sample equator attains H c \\ vor- 
tex sections penetrating from the top and the bottom 
ridge join, and the vortex line enters the bulk. A key ex- 
periment demonstrating the relevance of the geometrical 
barrier in Bi2Sr2CaCu20s+5 is the measurement of the 
penetration field of crystals with similar aspect ratio d/w, 
but different shape: a prism-shape d 15 i 31 or ellipsoidal^ 2 - 
crystal have a lower penetration field than a rectangular 
parallelepiped. Experiments reporting the accumulation 
of flux in the center of a superconductor sample imme- 
diately after first flux penetration cannot be taken as 
evidence for the presence of a geometrical barrier; such 
observations merely reflect the fact that the sample cen- 
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FIG. 1: (color online) Schematic diagram of vortex penetra- 
tion in (a) a "continuous" superconductor in which layeredness 
is irrelevant and (b) in a layered superconductor; both mate- 
rials are presumed to have a rectangular cross-section. Vortex 
lines (in blue) first penetrate at the corners, and then progress 
inwards to the sample equator. The direction of the screening 
current is indicated by the red circles. In case (a), there is no 
further barrier impeding the progression of the vortex lines 
(arrows); furthermore, the vortex currents run parallel to the 
lateral surface. Therefore, no effect of the Bean-Livingston 
barrier superimposed on that of the geometrical barrier is 
to be expected. In (b), the progression of the vortex line is 
mediated by the periodic entry of pancake vortices (blue rect- 
angles); each of these should overcome the BL barrier. The 
supercurrents associated with the pancakes run in a plane 
perpendicular to the lateral surface, therefore the geometrical 
barrier is enhanced by the Bean-Livingston barrier. 



ter corresponds to the location of minimal vortex free 
energy, a situation that stems from the macroscopic dis- 
tribution of the Meissner current and is therefore true 
regardless of the origin of the barrier. 

The flux density and current distributions accompa- 
nying vortex penetration into superconductors of rect- 
angular cross-section and varying thickness-to-width ra- 
tio d/w was numerically calculated by Brandt j 16 i 17 The 
penetration field obtained from these calculations is well 
described by the formula 

H p = H%>°° tanh (y/fid/w} . (1) 

The constant /3 was found to be ~ 0.36 for infinite strips 
(bars), and ~ 0.67 for rectangular cylinders; in the ab- 
sence of a BL barrier or vortex pinning, H^°° = H c \. 
Tentative agreement of Eq. (pQ) with H p of micron-sized 
single crystalline Bi2Sr2CaCu208+<5 was obtained for a 
few values of d/w and temperature T by Wang, et al. 
ji£ while Cabral et al. could fit their measured hystere- 
sis loops to those calculated in Ref. [13, but only in the 
temperature interval 50 - 60 KJ^ 

The purpose of the present paper is twofold. In a first 
step, we shall verify the dependence of the field of first 



flux penetration in square-shaped Bi2Sr2CaCu208+<5 sin- 
gle crystals on a wide range of aspect ratios d/w, over 
the entire temperature range from 20 K to T c . This is 
achieved by cutting the same crystal into successively 
smaller squares. These measurements establish the man- 
ner by which the Bean-Livingston barrier and magnetic 
relaxation influence H p , and yield the geometrical fac- 
tor f3 for a crystal of square shape. The results are then 
used to verify whether any effect of crystalline orienta- 
tion on Meissner screening and first flux penetration can 
be measured in single crystalline Bi2Sr2CaCu20s+5. 

The expected orientation dependence 24 should ac- 
tually be largest in layered superconductors such as 
Bi2Sr2CaCu208+5. In thin rectangular crystals of "con- 
tinuous" superconductors such as YBa2Cu307_5, there 
is no potential barrier stopping the progression of vor- 
tex lines into the interior once the geometrical barrier 
is overcome (Fig. [1]). Furthermore, the boundary con- 
dition imposes that the vortex current runs parallel to 
the crystal lateral surface. In this situation, any effect of 
the ABS quasi-particle current should be confined to the 
volume near the crystal ridges. In layered superconduc- 
tors, the progression of a vortex line is mediated by the 
successive penetration of pancake vortices. The entry of 
each pancake is counteracted by a potential barrier, since 
the pancake currents are always in a plane perpendicular 
to the lateral surface. These arguments hold for HTSC 
as long as the applied field is not exactly parallel to the 
(a, b) plane, since the high crytalline anisotropy impedes 
vortex penetration along the c-axis. A possible contri- 
bution of ABS to the Bean-Livingston barrier in layered 
HTSC will come not only from the crystal ridges, but 
from all fiat regions of the lateral surface that are ori- 
ented perpendicular to a nodal direction (a = 45°), and 
all wedge-shaped regions with bisector along the nodal 
direction^ In what follows, we find a tiny enhancement 
of Hp for T < 50 K, i.e. , in the regime in which flux pen- 
etration is delayed by the Bean-Livingston barrier. This 
shows that, if any enhancement of the Bean-Livingston 
barrier due to ABS exists, it is either much smaller than 
theoretically expected, or marred by physical damage at 
the crystal edges. At higher temperature, no enhance- 
ment is found. 



II. EXPERIMENTAL DETAILS 

A rectangular Bi2Sr2CaCu20s+5 single crystal is cut 
along its a and b axes from a seed rod grown by the trav- 
elling solvent floating zone method. 34 The crystal orien- 
tation was determined from the growth direction of the 
rod^ and verified by the precession X-ray diffraction 
method^ 

Flux penetration below 80 K was imaged using di- 
rect magneto-optical imaging (MOI)^ while for higher 
temperatures the Differential Magneto- Optical (DMO) 
metho d 37 ! 38 was used. In MOI, a ferrimagnetic garnet 
indicator with in-plane anisotropy is placed on top of the 
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FIG. 2: (a) DMO image of Meissner screening by the initial 
Bi2Sr2CaCu2 08+<5 single crystal, at 86.0 K and H a — 0; the 
field modulation AH a = 1 Oe. (b) Cutting scheme. The 
black rectangle represents the initial crystal, white numbered 
lines indicate successive cuts, in order. 



sample under study, and observed using a polarized light 
microscope. In all cases, the magnetic field, of magnitude 
H a , is applied parallel to the smallest crystal dimension 
and the c-axis. The presence of a non-zero perpendicular 
component B of the magnetic induction is revealed, by 
virtue of the Faraday effect of the garnet, as a non-zero 
intensity of reflected light when the polarizers of the mi- 
croscope are (nearly) crossed (see Fig. 2]). Images are ac- 
quired by an automated procedure, which also ramps the 
magnetic field, at the same rate, fiodH a /dt = 5 x 10 -4 
Ts _1 , for all experiments. In DMO, an image is taken 
at an applied field H a ; the field is then decreased to 
H a — AH a (with AH a = 1 Oe), and another image is 
acquired. The latter image is then subtracted from the 
former. The procedure is repeated 100 times, whence 
the 100 differential images are averaged (see Fig. [2^, and 
the inset to Fig. [5]). The use of the DMO technique was 
limited to temperatures above 82 K, at which the field 
modulation AH a is not screened once the magnetic field 
exceeds H p . The field of first flux penetration was deter- 
mined visually from both the direct and the differential 
magneto-optical images, taken at successive values of H a , 
from the appearance of a non-zero B within the sample 
contours (Fig. [4]). A more precise value could be obtained 
by calibrating the local intensity of the light; this yields a 
map of the local induction from which curves of the local 
magnetic hysteresis were deduced. Similar curves were 
obtained using a micron-sized Hall probe array* 1 ^ Exam- 
ples of these are shown in Fig. [6] The penetration field 
was then estimated from the peak in the hysteresis curves 
(jjiodH a /dt = 1 x 10 -4 Ts _1 for these measurement). 



In order to avoid any artefacts due to tiny vari- 
ations in sample composition 39 or doping, the same 
Bi2Sr2CaCu208+5 crystal is chosen for all further exper- 
iments. The selected crystal, with T c = 88.7±0.5K, had 
initial dimensions (I = 1400) x (w = 500) x(d= 28) /im 3 , 
and was checked using the differential magneto-optical 
(DMO) metho d 37 ! 38 to be as homogeneous as possible: 
macroscopic and mesoscopic inhomogeneities are known 
to affect H p i 20 ^ Nevertheless, some inhomogeneities 
arising from the crystal growth method remained, 41 Fluc- 
tuations in crystal growth rate result in small spatial 
variations of the composition, with a arc-shaped topology 
reminiscent of the floating-zone meniscus (inset of Fig.[5j). 
The associated spatial variation of the local T c leads to 
preferential flux penetration in the arc-shaped regions, 
above the temperatures marked by arrows in Fig. [5j 

For the measurements of H p as function of d/w, the 
same roughly square shape is retained to eliminate the 
influence of sample shape, or ambiguities related to the 
difference between length and width. Using a W wire saw 
with 1 /xm SiC abrasive grains, the initial crystal, with 
d/w = 0.045, is progressively cut into five ever smaller 
squares, down to d/w = 0.47, according to the scheme 
of Fig. [2}d. In the final step, one of the intermediate 
squares is cut at 45° with respect to the crystal axes; it 
had an aspect ratio of d/w = 0.12. Fig. [3] shows scan- 
ning electron micrographs of the edges obtained following 
the cuts. Typically, different layers of the crystal split 
at slightly different distances x from the crystal center, 
yielding an edge roughness ((x 2 ) — (x) 2 ) 1 / 2 ~ 2 — 3 /im. 
The edges of individual layers are smooth on the scale of 
the penetration depth, ~ 200 nm, with occasional protru- 
sions. The local orientation of the edge meanders about 
the global orientation of the crystal, with a rms devia- 
tion of 10-25°, depending on the investigated edge. No 
systematic morphological differences could be observed 
between [100] or [010] "anti-nodal" and [110] "nodal" in- 
terfaces. The crystal edges therefore seem to be suitably 
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FIG. 3: Scanning electron micrographs of crystal edges cut 
at 90° (a-c) and at 45° with respect to the principal [100] 
and [010] crystal axes (d-f). Images (a,b;d,e) have the same 
magnification, as do images (c,f). 
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FIG. 4: (a) Flux penetration into the crystal with d/w = 0.33, 
delimited by lines "3", "3", and "2" in Fig.^b), at 70.0 K, and 
two fields below (fioH a — 4.9 mT) and above (/JoHa — 7.0 
mT) H p , respectively, (b) Scanning electron micrograph of 
the crystal cut at 45°, with arrows showing the areas en- 
larged in Fig. [3] (d,e); the subsequent three images show flux 
penetration into this crystal as visualized by MOI at 80.0 K 
and H a = 10.5 Oe < H p , 15.0 Oe « H p , as well as H a = 21.9 
Oe> H p . 



described by the model advocated in Ref. l33l . 
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FIG. 5: (color online) Enlarged view of the penetration field 
data of Fig. [3 in the temperature regime close to T c . Closed 
symbols denote data obtained by Hall probe magnetome- 
try, open symbols are data measured using magneto-optical 
imaging. Data above 82 K are obtained using the differen- 
tial magneto-optical method. Arrows show the temperature 
above which flux penetration is influenced by crystal inhomo- 
geneity. For example, the inset shows flux penetration in the 
initial crystal of Fig. [2] (d/w = 0.045) at 87.0 K, measured 
by the DMO method at a field of 3 Oe, and a field modula- 
tion AH a — 1 Oe. The observed heterogeneity leads to the 
"collapse" of H p at T = 86.8 K. Lines are guides to the eye. 



III. RESULTS 

The main panel of Fig. [5] shows H p of the different 
crystals cut from the initial sample in the vicinity of T c . 
One systematically observes a linear decrease of H p (T) 
upon approaching T c , followed by a steep descent to H p = 
(see also Ref. 42). DMO imaging shows that this drop is 
the consequence of flux penetration into areas of lower T c , 
and is therefore the hallmark of sample inhomogeneity. A 
similar depression of the field of first flux penetration was 
obtained by introducing heavy ion-irradiated regions. 40 
In what follows, we shall discard data taken in the high 
temperature regime of inhomogeneous flux penetration. 

Fig. collects the penetration fields for all tempera- 
tures and all thickness/ width ratios. As expected from 
the demagnetizing effect, the field of first flux entry 
strongly depends on the aspect ratio of the sample. Wide 
thin platelets have a lower penetration field than narrow 
thick ones. This shows that the local induction near the 
edge is mainly determined by the shape of the crystals, 
and that edge roughness and protrusions (Fig. [3]) do not 
play a determining role in defining H p . The measured 
dependence of H p on d/w, presented in Fig. [8j is very 
well described by Eq. (pQ). However, the experimental 
values f3, ranging between ~ 2 — 6 (inset to Fig. [7j), 
are significantly larger than expected. The prefactor 
Hp ,oc determined from the fits is presented in Fig. [9j 



Above T ~ 45 K, it corresponds to H c i(T) of the d- 
wave superconductor— in the limit d/w — »• oo, with the 
correct^ zero-temperature extrapolated H^O) = 180 
Oe and T c = 88.7 K as obtained from DMO. Previ- 
ously attempted comparisons to the two-fluid relation^ 



H$> oo (T) = H$> oo (0) 



yield an anomalously 
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low value H cl (0) « 100 Oe. 

We note the observation of thermally activated flux 
entry at all temperatures, indicative of a role of the BL 
barrier. In particular, the width of the hysteresis loops 
of Fig. [6] progressively narrows as T increases, until the 
loops are all but closed for T < T c . As in Ref. l8|, the 
temperature dependence of the penetration field in the 
regime 40 - 50 K is in qualitative agreement with the 
theory of Ref. 0, which has 
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The relaxation rate is observed to be S ~ —0.03. How- 
ever, if one attempts a quantitative fit by inserting this 
value, which should correspond to S = dH p /dlnt = 
UbT /eq(T)s, into Eq. ([2]), one obtains i^ p -values that, 
above 50 K, fall vastly short of the measured ones. There- 
fore, vortex entry above this temperature is not limited 
by the BL barrier (see also Ref. M). Below T - 30 K, 
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the shape of the hysteresis loops and the observation of 
remnant flux clearly indicate the effect of flux pinning. 
H p shows an associated steep exponential increase as T 
is lowered. 

The fits to Eq. (pQ) yield a set of master curves allowing 
us to assess the dependence of H p on the relative orien- 
tation of the sample edges with respect to the principal 
crystal axes. For this, we present, in Fig. [9j the H^ 5 
data measured on the sample with d/w = 0.12, cut at 
45°. These are compared to the interpolated values of 
H® for a hypothetical sample with the same aspect ratio, 
cut parallel to the principal axes (a = 0). At temper- 
atures above 50 K where the Bean-Livingston barrier is 
inoperative^ we find H^ 5 to be slightly smaller than H®. 
However, in the regime where the Bean-Livingston bar- 
rier impedes pancake vortex entry into the crystal, there 
is a small enhancement H^ 5 / H® of up to about 1.25. This 
enhancement increases with decreasing temperature. 



IV. DISCUSSION 

We first discuss the enhancement of the penetration 
field by the BL barrier and the effect of sample shape 
for a = 0. The dependence of H p on d/w is well de- 
scribed by Eq. (pQ) at all temperatures, not only at those 
where the barrier is of purely geometric origin. This is a 
consequence of the macroscopic distribution of the Meiss- 
ner current, which depends only on sample shape and the 
value of the penetration depth. At temperatures below 50 
K, pancake vortex entry is hindered additionally by the 
BL barrier, which leads to a larger penetration fields 
This enhancement is manifest through the increase of 
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FIG. 6: (color online) Hysteresis loops of the magnetic flux 
density in the center of the crystal with d/w — 0.052, mea- 
sured by a microscopic Hall probe array. The inset shows an 
expanded view of the loops at high temperature. 



f3 as temperature is lowered. As temperature increases, 
pancake vortex penetration is accelerated due to ther- 
mal activation over the BL barrier. However, the field of 
first flux penetration can never drop below that of the 
geometrical barrier, since below this field H < H c \ at 
the crystal equator and the existence of vortices in the 
superconductor is thermodynamically unfavorable. Once 
Hp is overcome, thermal activation of pancakes leads to 
the rapid increase of the local induction at the sample 
boundary to values above H&, and the concommitant 
closing of the magnetic hysteresis loops. 

Even in the high temperature regime where the BL 
barrier has no effect whatsoever, the experimentally 
found /3 is much larger than that found in the calculations 
of Refs. 16 and [I?]. Such a large value of f3 expresses the 
fact that in thinner crystals, the measured penetration 
field is larger than what would be expected from calcula- 
tions including the effect of the geometrical barrier only. 
An explanation is that, apart from opposing pancake en- 
try from the side surfaces, the Bean-Livingston barrier 
opposes initial flux penetration across the top and bot- 
tom sample ridges, a finite size-effect that is expected to 
be more important for thinner crystals. Furthermore, j3 
is enhanced because of the square shape of the sample 
surface perpendicular to the field direction. Although 
this is manifestly closer to the disk than to the strip ge- 
ometry, the corners of a square superconductor are re- 
sponsible for enhanced screening with respect to the disk, 
yielding a larger penetration field. The effect is similar 
to that of the ridges of a platelet of rectangular cross- 
section, which are responsible for additional screening 
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FIG. 7: (color online) The penetration field as evaluated 
from MO-imaging (open symbols) and Hall probe mea- 
surements (full symbols). The aspect ratio is d/w — 
0.47, 0.33, 0.12(45°), 0.088, 0.052, 0.045, from the upper to 
the lower curves. The inset shows the experimentally deter- 
mined factor p as function of temperature. 
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FIG. 8: (color online) Double- logarithmic plot of the penetra- 
tion field of single crystalline Bi2Sr2CaCu2 0s as a function of 
aspect ratio for different temperatures. Open symbols show 
for the case where the edges are parallel to the principal 
axes. Drawn lines are fits to Eq. (pQ). Full symbols indicate 
Hp 5 (edges at a = 45° with respect to the principal crystal 
axes) . 



when compared to an ellipsoidal superconductor of the 
same d/w. By extrapolating to T c , at which all effects 
of the Bean-Livingston barrier have vanished, we find, 
experimentally, that j3 ~ 1 for a square superconduc- 
tor with rectangular cross-section. Note that taking an 
incorrect value of /? may result in large error when de- 
termining of the first critical field H c \ from H p of thin 
superconducting platelets. 

We now turn to the relative enhancement H^ 5 / H® 
of the field of first flux penetration for a = 45°, i.e. 
where the crystal edges are perpendicular to the nodal 
directions of the superconducting order parameter in 
Bi2Sr2CaCu208+(5. The observation of an enhancement 
in the sole temperature regime where the BL barrier is 
relevant is in qualitative agreement with the theory of 
Iniotakis, Dahm, and Schopohl. The theoretically pre- 
dicted enhancement H^ 5 / H® ~ 1 + c//£ increases as 
temperature decreases (see Inset to Fig. [9]). Its physical 
origin is that, in superconductors where the order pa- 
rameter A undergoes a sign change on any given Fermi 
surface sheet, interference between quasiparticles with 
a relative phase difference of tt will result in the oc- 
curence of Andreev bound states (ABS) on those parts 
of the superconductor- vacuum interface that are perpen- 
dicular, or nearly perpendicular, to the /c-space direc- 
tion for which A vanishes. For Bi2Sr2CaCu208+5 with 
d x 2_ y 2-w&ve symmetry of the order parameter^ these 
"nodal" directions correspond to (±7r, ±7r), i.e. at a = 45° 
with respect to the principal crystalline axes. 26 The ABS 
have a spatial extent of several times £ in the direc- 
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FIG. 9: (color online) Plot of the penetration field Hp 5 of a 
Bi2Sr2CaCu2 08 crystal with the edges cut at a = 45° with 
respect to the main crystal axes ), compared to the inter- 
polated values Hp for a crystal with the same d/w = 0.12, 
but a = (q ). Also shown is Hp°°(T) for a crystal with 
d/w ^ oc, a = 0, extracted from the fits in Fig. [8] ( ); the 
drawn line is a fit to the theoretical temperature dependence 
of Hd for a d-wave superconductor,— with H c i(0) = 180 Oe 
and T c = 88.7 K. The inset shows the experimental (lefthand 
axis) and the theoretically expected enhancement (righthand 
axis) of the field of first flux penetration for a = 45° with 
respect to a = 0. 



tion perpendicular to the interface, and are accompa- 
nied by a local reduction of |A|. Translational invari- 
ance along the interface allows these states to carry cur- 
rent. When a magnetic field is applied to the super- 
conductor, this "anomalous" paramagnetic current adds 
to the diamagnetic screening current, and to the ordi- 
nary paramagnetic current resulting from the Doppler 
shift of quasiparticle states. 44 In the region near a nodal 
surface, the total Meissner current is therefore reduced 
with respect to surfaces that are perpendicular to the 
/c-space direction where the gap is maximum^ hence, 
the magnetic field should penetrate somewhat further 
into the superconductor^ The presence of a vortex 
near the nodal interface should also affect the Andreev 
bound states j 4 - 6 with the result of somewhat reducing the 
anomalous current. 24 The BL barrier on the nodal surface 
must therefore be evaluated by the addition of the stan- 
dard free energy contributions from the vortex- Meissner 
current repulsion and the attraction of the vortex to the 
surface; and the "anomalous" contributions arising from 
the extra repulsion due to the augmented field near the 
surface, and from the attraction between the vortex cur- 
rent and the bound state current j 2 - 4 The net result is an 
enhancement of the barrier maximum, and a shift of its 
location towards the bulk of the superconductor. The en- 
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hancement of the BL barrier should be maximum when 
the orientation of the superconductor-vacuum interface 
is perpendicular to the nodal direction (and nil when the 
interface is parallel to the nodal direction), and be tem- 
perature dependent. For example, at T/T c = 0.2, an 
increase of H p by a factor of 4 to 5 is expected on sur- 
faces oriented at a = 45°, with respect to the situation 
where a = 0.— 

The observation of an enhancement of H p for the crys- 
tal with the edges cut perpendicular to the nodal di- 
rections seems to confirm the prediction by Iniotakis et 
al. of an orientation-dependent H p due to a more im- 
portant effect of Andreev bound states in this case i 24 i 33 
The inset to Fig. [9] shows that the magnitude of the ex- 
perimentally measured effect is, however, far less than 
predicted. 24 The prime factor that may influence the ex- 
perimental result is the roughness of the edges cut by 
the wire saw. The undulating character of the crystal 
edges will entail a reduction of the local density of bound 
states, and should strongly diminish the total anomalous 
(paramagnetic) Meissner current. As for the observed 
wedge-shaped protrusions, recent work has shown even 
if the occurence of ABS depends sensitively on quasi- 
particle reflection at the boundary as well as the shape 
and precise orientation of the latter^ they are still ex- 
pected as long as the wedge bisector is aligned close to 
a nodal direction^ 3 - However, the localised quasiparticles 
near to such protrusions will not contribute to the total 
Meissner current, and will have no effect on the Bean- 
Livingston barrier. In summary, the area on which the 
Bean-Livingston barrier is effectively enhanced is much 
smaller than the total area of the lateral edges. Only 
a fraction of pancakes composing an entering vortex line 
will effectively be subjected to the influence of local ABS, 
entailing the decrease of H p 

Finally, even in those parts of the surface where ABS 
would be fully developped, pancake vortices enter the 
crystal by thermal activation over the barrier. The im- 
plies that pancakes are nucleated at a distance that 
is farther removed from the edge than the predicted 
iH c / '\[2H a ^ In Ref. 0, the magnitude and temperature 
dependence of the enhancement factor 1 + c//£ strongly 
depends on the distance of the vortex from the crystal 
edge. Any enhancement of this distance due to thermal 
activation would imply a much smaller modification of 
the Meissner current, and a significant reduction of the 



enhancement factor. 



V. CONCLUSIONS 

Summarizing, we have measured the field of first 
flux entry, H pi for samples of different aspect ratio 
and orientation, cut from the same optimally doped 
Bi2Sr2CaCu208 single crystal. We find that at all tem- 
peratures, the penetration field of the square-shaped 
crystals is satisfactorily described by Eq. (pQ). The nu- 
merical factor /3 is determined jointly by the shape of 
the sample surface perpendicular to the field direction 
and by the effect of the Bean-Livingston barrier, which, 
in layered superconductors, counteracts the penetration 
of individual pancake vortices. For the case of a square 
superconductor of rectangular cross-section, without the 
influence of a Bean-Livingston barrier, it is experimen- 
tally found that f3 ~ 1. 

Moreover, we have investigated the enhancement of H p 
when the edges are perpendicular to the nodal directions 
of the order parameter, i.e. at 45° with respect to the 
crystalline a and b axes. If discernible, the measured en- 
hancement is tiny and occurs only in the temperature in- 
terval in which the Bean-Livingston barrier increases the 
field of first flux penetration over what would be expected 
from sole geometrical effects. The results are discussed in 
the framework of possible Andreev bound states present 
at the nodal surfaces of the d-wave superconductor. 24 
Notably, the small value of the enhancement in compar- 
ison to theoretical prediction is ascribed to edge rough- 
ness and thermal activation of pancake vortices over the 
Bean-Livingston barrier. 
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